Mon. Not. R. Astron. Soc. 000,[T]|5](2006) Printed 5 February 2008 (MN IATpX style file v2.2) 



A near-infrared spectroscopic detection of the brown dwarf in the 
post common envelope binary WD 0137-349 
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ABSTRACT 



We present a near-infra red spectrum o f the close, detached white dwarf + brown dwarf 
binary WD 0137-349 jMaxted et alJl2006l) . that directly reveals the substellar companion 
through an excess of flux longwards of w 1.95/um. We best match the data with a white 
dwarf + L8 composite model. For ages ~ 1 Gyr, the spe c tral type of the cool secondary is 
in agreement with the mass determined by Maxted et al. (2006) from radial velocity mea- 
surements (0.053 ± 0.006Mq), and supports an evolutionary scenario in which the brown 
dwarf survived a previous phase of common envelope evolution which resulted in the forma- 
tion of this close binary. The brown dwarf is the lowest mass companion to a white dwarf 
yet found, and the lowest mass object known to have survived a common envelope phase. At 
1300 < T eS < 1400 K WD 0137-349B is also the coolest known companion to a white dwarf. 
At a separation a = 0.65Rq the hemisphere of the brown dwarf facing the 16, 500K white 
dwarf is being heated through irradiation. We discuss the possible effects of this additional 
heating, with particular relevance to those other close binaries with substellar companions, 
the hot Jupiters. We propose future observations to investigate the likely temperature differ- 
ences between the "day" and "night" sides of the brown dwarf. 
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1 INTRODUCTION 

Detac hed brown dwarf companion s to white dwarfs are rare 
I Farihi. Becklin & Zuckerman 2005). Proper motion surveys and 
searches for infra-red (IR) excesses have s o far f ound only three 
confir med examples: GD 165 (DA+L4, Becklin & Zuckerman 
I l988t). GD1400 (DA+L6/7, iFarihi & Christopher! 12001 
bobbie et al] 120051) . and WD 0137-349 JMaxted etaU 120061) . 
the subject of this paper. GD 165 is a widely separated system 
(120 AU); the separation of the components in GD 1400 is 
currently unknown. In contrast, WD 0137-349 is a close binary 
(orbital period P ~ 116 minutes). 

Optical spectra of the H-rich DA white dwarf WD 0137-349 
show a narrow H a emission line due to irradiation of the surface 
of the unseen companion. Radial velocities measured from this 
line and the white dwarf's intrinsic Ha absorption line allowed 
lMaxtedetail f2006) to determine the mass ratio of the system. Us- 
ing the white dwarf mass (0. 39 ± 0.035Mp) ), derived from an anal- 
ysis of its optical spectrum, Maxted et al. 1 2006 ) then determined 
the mass of the companion to be 0.053 ± O.OOBM©. This is well 
below the limit of 0.075Mq commonly used to distinguish stars 
from brown dwarfs. The substellar nature of WD0137-349B was 
reinforced by an analysis of its 2MASS near-IR JHK fluxes (Fig- 



ure 2, lMaxtedetall2 006). The J and H band photometry can be 
fit by a model consistent with the white dwarf alone. There is a 
slight excess of flux at Ks, which can be matched by a model con- 
sisting of the white dwarf plus a brown dwarf spectral type mid-L 
or later. This is consistent with the radial velocity determined mass 
measurement for WD0137-349B. 

Therefore, WD 0137-349 is the first close, detached binary to 
be discovered containing a confirmed substellar companion. The 
brown dwarf must have survived a previous phase of common en- 
velope (CE) evolution during whi ch it was e ngulfed by the red gi- 
ant progenitor of the white dwarf iPolitano 2004). The drag on the 
brown dwarf caused it to spiral in towards the red giant core from 
an originally much wider orbit. Some fraction of the orbital energy 
was released and deposited in the envelope, which was ejected from 
the system, leaving a close binary. Simple physical arguments sug- 
gest that low mass companions less than some limit m C rit will be 
evaporated during the CE phase. WD0137-349B places an upper 
limit on m cr it. 

Alternatively, WD0137-349B may have originally been a 
planet which accreted a substantial fraction of its mass during the 
CE phase <Livio & Soke3ll983h . In this scenario, WD0137-349B 
would be expected to have an effective temperature T e « > 2000K, 
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and a spectral type (early L) commensurate with the cooling age of 
the white dwarf (~ 0.25 Gvrl. lMaxted et alj |20o3) concluded that 
this model was probably not applicable to WD 0137-349 since the 
2MASS fluxes are inconsistent with a companion earlier than mid- 
L. However, they suggested that a definitive test would be provided 
by IR spectroscopy to directly determine an accurate spectral type 
for the brown dwarf. 

In this paper we present a near-IR spectrum of WD 0137-349. 
In Section 2 we discuss the observations and data reduction, in Sec- 
tion 3 we present our analysis of these data, and in Section 4 we 
discuss the results and comment on the spectral type of WD 0137- 
349B. 



2 OBSERVATIONS AND DATA REDUCTION 

We were awarded Director's Discretionary Time at the 8m Gemini 
South telescope in November 2005, to obtain a near-IR spectrum 
of WD 0137-349 with the Gemini Near-IR Spectrometer (GNIRS, 
lEliasetalll 998) in programme GS-2005B-DD-9. We selected the 
cross-dispersed mode, using the short camera, the 32 lines/mm 
grating centred at 1.65pm, and the 0.3" (2 pixels) slit, which gives 
a resolution R — 1700. In this observing mode, the entire near-IR 
region from ~ 0.9/im to ~ 2.5/im is covered in a single obser- 
vation with excellent transmission across almost the whole wave- 
length range. There is no inter-order contamination as in long slit 
mode over a single atmospheric window, and this mode is much 
more efficient than executing separate H and K band observations. 

The observations were conducted in service mode on the night 
of 2005 November 22nd. Our observing condition requirements 
were met: 70%-ile image quality (seeing < 0.6" at J) and 50%-ile 
cloud cover (i.e., photometric). The observations were obtained at 
a mean airmass of 1.12. We took 20 x 120s exposures, using two 
"nod" positions along the 6" long slit giving a total exposure time 
of 40 minutes. This did not cause any problems of order overlap 
during the data reduction, since the target is a point source and the 
seeing was good. With overheads for detector readout etc., the ob- 
servations lasted for ~ 55 minutes, covering « 47% of the binary 
orbit. An Al V telluric standard was also observed (HD 10538) 
both before and after the target, at a similar airmass. 

The data were initially reduced using the GNIRS sub-package 
of the Gemini IRAF package. Briefly, a correction was first applied 
to the raw science, standard star and arc lamp spectral images for 
the s-distortion in the orders. The data were then flat-fielded, tak- 
ing care to flat-field each order with the corresponding correctly 
exposed flat. Subsequently, difference pairs were assembled from 
the science and standard star images and any significant remaining 
sky background removed by subtracting linear functions, fitted in 
the spatial direction, from the data. The spectral orders of the white 
dwarfs and the standard stars were then extracted and assigned the 
wavelength solution derived from the relevant arc spectrum. 

Unfortunately, the GNIRS IRAF routines proved inadequate 
for removing the telluric features from the target spectra and pro- 
viding an approximate flux calibration. Instead, at this stage we 
utilised standard techniques offered by software routines in the 
STARLINK packages KAPPA and FIGARO. Any features intrin- 
sic to the energy distribution of the standard star were identified by 
refere nce to a nea r-IR spectral a tlas of f undamental MK standards 
(e.g. IWallace & Hinklel Il997t IWallace et ail l200Ct iMever et alJ 
1998) and were removed by linearly interpolating over them. The 
spectrum of the white dwarf was then co-aligned with the spec- 
trum of the standard star by cross-correlating the telluric features 



present in the data. The science spectral orders were then divided 
by the corresponding standard star spectral orders and multiplied 
by a blackbody with the standard star's T c g , taking into account 
the differences in exposure times. 

We found that the flux levels of the three brightest spectral 
orders (number 3, covering ~ 1.95/im to ~ 2.5/im, number 4 cov- 
ering ~ 1.4/im to ~ 1.8/im, and number 5 covering ~ 1.2/xm 
to ~ 1.35/im) were well matched to each other. However, it was 
necessary to scale these fluxes by a single, constant normalisa- 
tion factor to obtain the best possible agreement between the spec- 
tral data and the J, H and it's photometric fluxes derived from 
the 2M ASS All Sky Data Re lease Point Source Catalogue mag- 
nitudes ISkrutskie et alJI 19951) where zero magnitude fluxes were 
taken from Zombeck 1 1990). The reduced GNIRS spectrum and 
2MASS fluxes are shown in Figure 1. 



3 ANALYSIS 

We have analysed the GN IRS near-IR spectrum of WD 0137-349 
following the method of Dob bie et alJ l2005t) . We calculated a 
pure-hydrogen synthetic white dwarf spectrum for the tempera- 
ture (T e ff = 16, 500 ± 500 K) and surface gra vity (log g = 
7.49 ± 0.08) determined by iMaxted et all J200d) . We used re- 
cent versions of the plane-parallel, hydrostatic, non-local thermo- 
dynamic equilibrium (non-LTE) atmosphere and spectral synthesis 
codes TLUSTY (V200: lHubenv & LanzkOOll and SYNSPEC (v48; 
ftp:/tlusty.gsfc.nasa.gov/synsplib/synspec). The synthetic spectral 
fluxes have been normalised to the white dwarf's V magnitude 
(15.33 ± 0.02). The synthetic white dwarf spectrum is shown over- 
plotted the GNIRS data in Figure 1. 

Examination of Figure 1 reveals a clear difference in the over- 
all shape and level of the synthetic white dwarf spectrum and the 
observed data in the K band (below w 1.95/im). In contrast, we do 
not readily observe any spectral features typical of the energy dis- 
tributions of L or T dwarfs, e.g. Na I absorption at 2.20/im, CH4 or 
CO at 1.6/x and 2.3/^m respectively and H2O centred on 1.15, 1.4 
and 1.9/im. We estimate the signal-to-noise ratio at 2.2/im ps 12. 

To test whether the flux excess in the K band is due to the 
brown dwarf companion, and to constrain its spectral type, we 
have added empirical models for substellar objects to the white 
dwarf synthetic spectrum and compared these composites to the 
near-IR spectrum (Figure 1). The empirical models have been con- 
structed using the near-IR spectra of L and T dwarfs presented by 
iMcLean et alJ 1200 3f) . In brief, the data have been obtained with 
the NIRSPEC instrument on the Keck Telescope, cover the range 
0.95/^m — 2.31/im with a resolution of \/8\ ~ 2000 and have 
been flux calibrated using J, H and Ks photometric fluxes derived 
from the 2MASS magnitudes as described by Mc Lean et ailfeOOlh . 
To extend these data out to 2.4/im, our effective red limit, we have 
appended to them sections of U KIRT CGS 4 spectra of late-typ e 
dwarfs obtained bv lLeggett etailfeOOll) and lGeballe et alJ 12002). 

The fluxes of the empirical models have been scaled to a level 
appropriate to a location at d=10pc using the 2MASS J magni- 
tudeof each late-type object and the polynomial fits o flTinnev et alJ 
1 2003) to the Afj versus spectral type for L0-T8 field dwarfs/brown 
dwarfs. Subsequently, these fluxes have been re-calibrated to be 
consistent with the distance of WD 0137-349, as derived from its 
measured V magnitude, effective temperature and theoretical My 
and radius. 

In Figure 1 we compare the observed GNIRS near-IR spec- 
trum of WD 0137-349 to a range of these empirical composite 
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Figure 1. The GNIRS spectrum of WD 0137-349 from 1.2/im — 2.4/im compared to a synthetic white dwarf spectrum from a pure-H model atmosphere 
normalized using the observed V band magnitude (solid line). Also shown are the synthetic white dwarf spectrum combined with spectra of known brown 
dwarfs scaled to the appropriate distance as follows (dashed lines, top-to-bottom): L0, L6, L8 and T5. Note the telluric water bands between 1.35/xm — 1.42/xm 
and 1.8/jm — 1.95/im have been omitted. 
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white dwarf + brown dwarf models (WD+LO, L6, L8 and T5). 
We find the best match to the data is provided by a white dwarf + 
L8 composite model. The lack of an obvious CO edge at 2.3/im 
supports this conclusion. 



4 DISCUSSION 

We can estimate the temperature of the brown dwarf from its 
spectral type. From astrom etric and photometric observations of 
L dwarfs Vrba et al. (2004) give the mean effective temperature of 
an L8 dwarf as 1390K, although they caution that their derived val- 
ues should be treated as schematic results only, and note that late 
L and early T dwarfs appear to occupy a fairly narrow temperature 
range from 1200 < T e g < 1550 K. The temperatures computed by 
Golim owski et alj 120041) for a small sample of late L and early T 
dwarfs are in agreement with these estimates. For example, they 
estimate the effective temperature of the L8 dwarf Gl 584C, whose 
age is constrained by its main sequence companion to 1 — 2.5 Gyr 

iKirkpatrick et all200ll) . as 1300 < T e S < 1400 K. 

Using the evolutionary models of Burrows et alj 1 19971) and 
their on-line calculator 1 , we find that for a 0.053Mq brown dwarf 
to have a temperature in the range 1300 < T e g < 1400 K it must 
be ~ 1.0 Gyr old. The COND models of lBaraffe et ai"h2003l) also 
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suggest that this mass and temperature range are consistent only for 
a cooling age ~ 1.0 Gyr. 

The low temperature of WD 0137-349B therefore supports the 
conclusion of M axted et alj 120061) that it is an old brown dwarf 
that survived a phase of CE evolution some ~ 0.25 Gyr ago, and 
does not support the rather more exotic evolutionary scenario of 
iLivio & Sokeri n*983) in which a lower mass planet or brown dwarf 
accreted substantial mass during the CE. In that case the tempera- 
ture would be substantially higher (T e fi > 2000 K) and the spectral 
type would be early L. 

If the brown dwarf is ~ 1.0 Gyr old, then the main sequence 
lifetime of the white dwarf progenitor was » 0.75 Gyr, in which 
case the progenitor's mass was just over 2M0. During the CE 
phase, some fraction of the orbital energy released, qce, will be 
deposited as kinetic energy in the envelope, which is ejected from 
the binary. A progenitor mass of 2MQ requires a value of qce > 2 
to explain the formation of the WD 0137-349 system. However, we 
note that the cooling age of the brown dwarf does not necessarily 
set an upper limit on the system age, since it may have been re- 
heated during the immediate post-CE phase when the white dwarf 
was very hot. The system may be older, in which case the white 
dwarf progenitor mass and the required value of qce are lower. 

The brown dwarf hemisphere facing the 16, 500K white 
dwarf is currently intercepting w 1% of its light. This irra- 
diation is seen as weak Ha emission in the optical spectrum. 
Therefore, WD 0137-349 is an interesting system for studying 
the effects of irradiation on the atmosphere of a substellar ob- 
ject, and can potentially be used as a comparison for differ- 
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ent theoretical models of the effects of irradiation on lower- 
mass hot Jupiters. The heated atmospheres of the transiting plan- 
ets HD 209458b and TrES-1 have both been detected by the 
Spitzer Space Telescope in the m id-IR during secondary eclipse 
toeming et all2005llCharbonneau et alj|2005t) . Irradiation can in- 
crease the photospheric temperature by an order of magnitude 
compared to isolated planets, decrease the cooling rate and al- 
ter the planet's radius and atmospheric structure (e .g. lGuillot et alJ 
1 19961 [Baraffe et al] |2003l [Arras & Bildsterl EoOa) . Severe irradi- 
ation could even lead to atmospheric evaporation ( Baraffe et al. 
|2004) . for which evidence has been found thr ough the discovery 
of an extended atmosphere for HD 209458b l Vidal-Madiar et al. 
120031 1 Vidal-Madiar et al1l2004l) . In these synchronously rotating 
systems, substantial temperature differences are expected between 
the "day" and "night" sides, possibl y leading to strong winds 
transporting heat to the "night" si de IShowman & GuilToll2002l 
iBarman.Ha uschildt & Allarci 2005t). S imilar effects may occur fo r 
more massive brown dwarfs (Hubenv. Burrows & Sudarskv 2003). 

Neglecting the contribution from the intrinsic luminosity of 
the brown dwarf, and assuming it is in thermal equilibrium with the 
stellar radiation, we can use the relation given bv lArras & Bildster] 
(2006) to estimate the "equilibrium" temperature of the irradiated 
hemisphere: 

T ecl = T*(R*/2a) 1/2 (1) 

Taking the white dwarf radius (Rwd = 0.0186Rpi) and bi- 
nary separation (a = O.65R0) given by iMaxted et alJ 120061) . 
we estimate the equilibrium temperature of the heated face of the 
brown dwarf to be T cq w 2000K. For ages ~ 1 Gyr, the spectral 
type of an object of this temperature would be early L. Clearly, 
we are not seeing a hemisphere at this temperature in our data. 
The GNIRS spectrum was obtained across 47% of the orbital pe- 
riod, so it is slightly possible that the heated face has not been ob- 
served. However, the 2MASS photometry are also not in agree- 
ment with a temperature this high, although the phase at which 
they were obtained is unknown. Eq. (1) assumes zero reflection of 
the stellar photons, and should be multiplied by (1 — A) 1 ^ 4 for 
an albedo A (assuming isotropic emission from the brown dwarf). 
Substituting our estimated temperature of WD 0137-349B from the 
GNIRS spectrum (T cB = 1300K) for T ea in Eq. (1), we find a high 
albedo A ~ 0.8. For comparison, Charbo nneau et al] 120051) find 
that TrES-1 appears to be absorbing a high fraction of the incident 
stellar flux (A ~ 0.3), although we note that the spectrum of the 
light falling on the surface of WD0137-349B is very different to 
that falling on a hot Jupiter. 

None-the-less, there is likely some heating of the hemisphere 
facing the white dwarf. A fraction of the observed flux could orig- 
inate from the heated face. This could mean that the underlying 
brown dwarf temperature might be lower than 1300K, and the de- 
duced age should then be regarded as a lower limit. 

We suggest that time series K-band photometry should reveal 
variability due to temperature differences between the "day" and 
"night" sides. Spitzer mid-IR observations should also reveal such 
temperature variations, allowing us to characterize the photomet- 
ric properties of the brown dwarf WD 0137-349B across the whole 
orbit. 



5 CONCLUSIONS 

We have obtained a near-IR spectrum of the close, detached white 
dwarf + brown dwarf binary WD 0137-349 that reveals the substel- 



lar companion through an excess of flux longwards of ~ 1.95/^m. 
We best match the data with a white dwarf + L8 composite model. 
This is the lowest mass, coolest and latest spectral type compan- 
ion to a white dwarf yet discovered. For a cooling age ~ 1.0 Gyr, 
the temperature of the secondary (130 < T e s < 1 400 K) is in 
agreement with the mass determined by M axted etalJfeOOd) from 
radial velocity measurements (0.053 ± 0.006Mq), and supports an 
evolutionary scenario in which the brown dwarf survived a previ- 
ous phase of CE evolution ~ 0.25 Gyr ago, which resulted in the 
formation of this close binary. Th e exotic alternative evolutionary 
scenario of iLivio & SokeJ 1 1983l> . in which a lower mass planet 
or brown dwarf accreted substantial mass during the CE phase, is 
not supported. If the brown dwarf is ~ 1.0 Gyr old, then the main 
sequence lifetime of the white dwarf progenitor was « 0.75 Myr 
and its mass was just over 2M©, requiring a value of the CE pa- 
rameter QfCE > 2 to explain the formation of the system. How- 
ever, the brown dwarf may be older if it was re-heated by the very 
hot white dwarf immediately after the CE phase. The brown dwarf 
hemisphere facing the 16, 500K white dwarf is currently intercept- 
ing « 1% of its light, but does not appear to be substantially heated, 
although our spectroscopic observations only cover « 0.5 of the 
orbital period. 
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